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Summary 


The  objective  of  this  research  project  is  to  investigate  the  anomalous  normal  incidence 
absorption  observed  in  n-type  direct  gap  quantum  well  infrared  photodetectors  (QWIPs) 
fabricated  from  the  GaAs/AlGaAs  and  strained  layer  InGaAs/AlGaAs  materials  systems 
grown  on  GaAs  substrates  as  well  as  developing  multi-color  QWIPs  for  3-5  pm  MWIR 
and  8-14  pm  LWIR  applications.  In  this  research  project,  we  have  conducted  theoretical 
and  experimental  studies  on  normal  incidence  absorption  in  «-type  direct  gap  QWIPs 
formed  on  AlGaAs/GaAs  and  InGaAs/AlGaAs  materials  systems.  A  theoretical  model  for 
calculating  the  TE  and  TM  transition  ratios  for  the  symmetrical  LWIR-  and  step  MWIR- 
QWIPs  has  been  developed  in  this  work.  The  role  of  strain  on  the  normal  incidence 
absorption  has  been  investigated.  We  have  also  calculated  the  absorption  coefficients  for 
several  n-type  QWIPs,  including  both  the  high  -strain  MWIR  QWIP  and  LWIR  TC- 
QWIPs.  Analysis  of  simplified  QWIPs  for  low  photon  background  applications  has  been 
made  in  this  work.  We  have  also  studied  the  quantum  confined  Stark  effect  in  a  weakly 
coupled  double  quantum  wells  structure.  Experimental  studies  on  three  QWIP  samples 
with  different  device  structures,  processing  parameters  and  experimental  conditions  have 
been  made  to  investigate  the  various  effects  on  the  normal  incidence  absorption.  Several 
new  simplified  QWIPs  have  been  designed  and  processed  which  show  strong  evidence  of 
normal  incidence  absorption.  Finally,  a  novel  four-stack,  four-color  QWIP  has  been 
designed,  fabricated  and  characterized.  Several  journal  and  conference  papers  have  been 
published  and  presented  from  this  work. 
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1.  Theoretical  Study  of  Normal  Incidence  Absorption  in  «-type  Direct  Gap  QWIPs 


1.1  Introduction 

Ever  since  West  and  Eglash  observed  experimentally  a  large  oscillator  strength  of 
the  intersubband  transition  in  a  series  of  multiquantum  wells  [1]  and  Levine  et  al. 
demonstrated  the  first  quantum  well  infrared  photodetector  (QWIP)  using  GaAs/Al^Ga,. 
^As  multiquantum  well  structure  [2],  a  considerable  amount  of  works  have  been  reported 
in  this  field  [3,4»5].  Taking  advantages  of  band  gap  engineering  and  the  maturity  of  the 
III-V  semiconductor  materials  growth  and  processing  technology,  large  uniform 
GaAs/AlGaAs  focal  plane  arrays  (FPAs)  with  low  noise  equivalent  temperature 
difference,  high  resolution,  high  production  yield,  high  speed,  low  cost  and  low  power 
consumption  have  been  demonstrated  in  recent  years. 

In  the  development  of  QWIPs,  an  important  issue  is  normal  incidence  absorption. 
It  is  well  known  that  normal  incidence  absorption  is  forbidden  in  n-type  direct  gap 
QWIPs  due  to  quantum  mechanical  selection  rules.  As  a  result,  it  is  a  common  practice  to 
use  Brewster  angle  or  the  45°  bevel  angle  incidence  configuration,  or  optical  grating  to 
couple  the  normal  incidence  light.  These  coupling  schemes  complicate  the  fabrication  of 
focal  plane  arrays  and  also  cause  cross  talk  between  the  adjacent  pixels.  Due  to  this 
reason,  there  has  been  a  continuous  interest  in  the  study  of  normal  incidence  absorption  in 
QWIPs.  P-type  QWIPs  with  normal  incidence  absorption  have  been  demonstrated  [6,7,8] 
and  the  mechanism  is  due  to  the  band  mixing  between  the  heavy-  and  light-  hole  states. 
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However,  due  to  the  large  heavy  effective  mass  and  low  mobility  of  holes,  the 
responsivity  of  jo-type  QWIPs  are  substantially  smaller  than  that  of  n-type  QWIPs. 
Normal  incidence  absorption  n-type  QWIPs  due  to  anisotropic  effective-mass  tensors 
grown  other  than  (100)  substrate  have  been  suggested  and  demonstrated  [9,10,11,12]. 
The  development  in  this  type  of  normal  incidence  absorption  QWIPs  is  limited  by  the 
difficulties  of  materials  growth  and  the  tilted  ellipsoidal  band  structure  also  causing  some 
complication  in  design  and  controllability. 

Normal  incidence  absorption  in  n-type  direct  gap  QWIPs  was  first  observed  by 
Peng  et  al.  [13]  in  InGaAs/InAlAs  system.  In  addition,  normal  incidence  absorption  in 
other  n-type  direct  gap  QWIPs  were  also  repeated  [14-21],  including  InxGa,.,(As/InyAl|. 
yAs  [13,14,16,17],  In,Ga,.,As/GaAs  [15,18],  GaAs/Al,Ga,.,As  [19,20,21]  and  In^Ga,. 
,(As/GaAs/AlyGa,.yAs  [20,21].  Theoretically,  several  mechanisms  have  been  proposed  for 
the  possibility  of  normal  incidence  absorption  in  n-type  direct  gap  QWIPs.  However, 
none  of  these  theories  is  able  to  explain  the  experimentally  observed  value. 

1.2  Summary  of  existing  theories  on  normal  incidence  absorption  mechanisms 

In  1981  Zalxxzny  pointed  out  that  in  narrow  gap  semiconductors  spin-orbit 
coupling  can  induce  spin-flip  excitations  by  IR  radiation  polarized  parallel  to  the  layer 
surfaces  [22].  Shik  considered  the  four-  band  Kane  model  and  found  that  both  spin-orbit 
coupling  and  the  variation  of  electron  effective  masses  in  space  can  induce  normal 
incidence  absorption  in  n-type  direct  band  gap  quantum  well  structure  [23,24].  Eight- 
band  Kane  model  has  been  considered  in  the  selection  rules  of  intersubband  transitions  in 
the  conduction  band  quantum  wells  by  Yang  et  al.  [25].  Yang  investigated  the  effect  of 
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the  spatial  variation  of  effective  masses  and  proposed  several  structures  to  enhance  the 
normal  incidence  absorption  in  «-type  direct  gap  QWIPs  [26, 27].  Peng  et  al.  developed  a 
generalized  14  -  band  model  to  discuss  normal  incidence  absorption  in  «-type  direct  gap 
QWIPs  [28].  They  showed  that  by  taking  into  account  the  coupling  between  r,®->r,5'', 
and  r,5'->r,5''  bands,  a  nearly  equal  TE  and  TM  absorption  can  be  obtained. 
However,  this  has  been  shown  to  be  inaccurate  lately  [29],  and  it  is  revealed  recently  that 
multiband  coupling  effect  is  not  sufficient  to  explain  the  observed  normal  incidence 
absorption  [30].  The  above  cited  studies  will  be  briefly  summarized  as  follows. 

1 .2. 1  One-  band  model  considering  spatial  variation  of  effective  mass 

Within  one-band  model,  if  the  effective  masses  of  electrons  in  the  quantum  well 
and  barrier  layers  are  taken  to  be  the  same,  the  in-plane  polarization  absorption  is 
prohibited.  This  is  due  to  the  fact  that  the  oscillator  strength  of  the  in-plane  (for  example, 
the  x-direction)  polarization  absorption  is  proportional  to  the  overlap  of  the  wave 
functions  of  initial  and  final  states;  since  different  states  are  orthogonal,  the  oscillator 
strength  will  be  zero,  i.e.. 


where  fx  is  the  x-polarization  absorption  oscillator  strength,  m*  is  the  effective  mass  of 
electron,  <j)j  and  (|)f  are  the  initial  and  final  state  wavefunctions,  respectively.  This  is 
schematically  shown  in  figure  1.2.1  (a). 
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By  taking  into  account  the  spatial  variation  of  the  electron  effective  mass,  the 
factor  of  the  inverse  mass  cannot  be  taken  out  of  the  integral  in  equation  (1),  and  hence 
will  have  a  nonzero  value,  i.e.. 


^0 


(2) 


This  is  schematically  shown  in  figure  1.2.1(b).  The  above  calculation  shows  that  the 
oscillator  strength  of  TE  absorption  is  generally  less  than  1%  of  that  of  TM  absorption 


[27]. 

1 .2.2  Kane  Model 


Another  way  to  study  normal  incidence  absorption  is  to  start  from  a  four-band  or 
an  eight-band  Kane  model.  Within  Kane  model,  the  band  coupling  and  spin-orbit 
interaction  cause  the  spin  degenerate  states  to  form  two  other  degenerate  states  which 
include  both  the  spin-up  and  spin-down  states,  as  shown  in  figure  1.2.2.  Normal 
incidence  absorption  can  occur  between  states  having  the  same  dominant  states  and  states 
having  different  dominant  states.  The  oscillator  strength  ratio  of  the  normal  incidence 
absorption  between  the  same  spin  group  to  the  z-polarization  absorption  is  almost  the 
same  as  the  ratio  mentioned  in  section  1.2.1,  which  is  less  than  1%.  The  oscillator 
strength  ratio  R  of  the  normal  incidence  absorption  between  different  spin  groups  to  the 
z-polarization  absorption  is  proportional  to  [23] 


|2 


3 


A 

2  ^ 


(3) 
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where  Eg  is  the  band  gap,  and  A  is  the  spin-orbit  splitting  of  the  well  material, 
respectively.  For  GaAs,  with  Eg~1.5  eV  and  A~0.34  eV,  we  obtain  R~0.004.  For  InGaAs, 
Eg~l  eV,  A~0.37  eV  and  one  finds  R~0.01.  Thus,  normal  incidence  absorption  cannot  be 
explained  using  the  four-  or  eight-  band  Kane  model. 

1 .2.3  Multi-band  coupling  effects 

Peng  et  al.  [28]  pointed  out  the  important  influence  of  multi-band  coupling  effect 
on  intersubband  transitions.  They  claimed  that  coupling  of  higher  conduction  band  Tu' 
with  the  valence  band  T,,''  and  lowest  conduction  band  T,®  can  cause  normal  incidence 
absorption,  as  shown  in  figure  1.2.3,  and  they  also  claimed  that  normal  incidence 
absorption  can  be  as  strong  as  the  z-polarization  absorption.  However,  their  analysis 
turned  out  to  be  incorrect  and  a  more  elaborated  numerical  calculation  showed  that 
normal  incidence  absorption  is  usually  negligible  within  this  multiband  scheme  [30]. 

1.3  Theoretical  Study  of  Normal  Incidence  Absorption  in  «-type  QWIPs 

1.3.1  One-  band  model  taking  into  account  spatial  variations  of  both  the  effective  mass 
and  refractive  index 

We  have  performed  a  theoretical  study  of  normal  incidence  absorption  on  several 
n  -  type  QWIP  structures  which  showed  obvious  normal  incidence  absorption  based  on  a 
one-band  model  taking  into  accoimt  spatial  variations  of  both  the  effective  mass  and 
refractive  index.  The  one-  band  model  works  well  for  the  study  of  TM  absorption.  On  the 
other  hand,  we  think  that  if  one-band  model  cannot  explain  normal  incidence  absorption, 
the  multi-band  model  may  not  be  able  to  explain  the  observed  normal  incidence 
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absorption  in  n-type  QWIPs  either.  Since  one-band  model  already  includes  the  influences 
of  valence  bands  and  higher  conduction  bands  through  the  effective  mass  of  electrons. 

The  envelope  function  in  the  z-  (the  growth)  direction  (p(z)  satisfies  the  following 
equation 


where  p  is  the  momentum  operator  of  electron  and  m*(z)  is  the  space- variant  electron 
effective  mass,  V(z)  and  e  are  the  potential  and  the  energy  of  electron  in  the  z-  direction, 
respectively. 

We  consider  the  symmetric  and  step  MQW  structure  with  one  bound  state  in  the 
wells,  as  shown  schematically  in  figure  1.3.1.  The  coordinates  adopted  in  the  calculation 
are  shown  in  figure  1.3.2.  The  wave  function  of  the  bound  state  can  be  written  as 

^4iexp(/(ri4  z <-/,./ 2. 

9’i  =  '  |z4z,../2.  (5) 

^,exp(_x-,z),  z>L,I2. 

where  Ly^,  is  the  well  width,  ki=(2m^,*  Ej/h^y^,  mi,*  and 

are  the  electron  effective  masses  in  the  barrier  and  well,  respectively,  Vi,  is  the  barrier 
height,  E I  IS  the  energy  level  of  the  bound  state. 

From  the  boundary  condition  one  obtains 

^1  =  BiCos(yt|£^/2.)exp(«:i£.,/2.)  (6) 

The  wave  function  for  the  continuum  state  with  odd  parity  can  be  written  £is 

.42sin(A:„r- J,),  z<-L../2. 

=  '  52sin(^j,,z),  |z|</,„,/2.  (7) 

[  /42Sin(A:2»z  +  (?2)-  z>I.r/2. 
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The  boundary  conditions  yield 


\ 


ftibklw 


Kmwkiby 


(8a) 


and  S2  can  be  determined  as 


(^2  =  arcsin 


2 


kih  Lw 


(8b) 


Similarily,  one  can  obtain  the  wave  function  of  the  first  even  parity  continuum  state  as 


<Pr 


A,coikv,^-S,)>  z<-Z,»/2. 

\z\<L,/2. 

A,<^ikv,^-^s,Y  z  >/,,./ 2. 


(9) 


where 


B,=  Ai 


cos2[^^l  + 


mbk3w]  .  if kiwLw 
mlkib)  ""  V  2 


(10a) 


and 


(5,  =  arccod 


2 


(10b) 


In  Eqs.  (7-10),  kih=(2mi,*(Ei-Vb)/^)^^^,  ki-^^;=(2my^*Ei/^)l^^ .  The  energy  level 
of  the  bound  state  and  the  central  energy  of  the  continuum  states  can  be  determined  by 
the  transfer  matrix  method.  Since  the  barrier  is  much  wider  than  the  well,  we  can  assume 

A2==A3^(2/Ly,+Lb)-^/2 

The  photon-electron  interaction  Hamiltonian  can  be  formulated  as 


4c 


1 

\m 


*P-^P  * 


J 
m  J 


(11) 
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where  q  is  the  wave  vector  of  the  optical  field.  Aq  and  e  are  the  amplitude  and  direction  of 
the  vector  potential. 


Incorporating  the  Hamiltonian  in  Eq.(ll)  into  the  transition  matrix  between  the 
initial  and  final  states,  the  oscillator  strength  due  to  in-plane  polarization  (chosen  to  be 
the  X  -  direction)  transition  van  be  expressed  as 


f 

J  if 


niw 


ha 


^  rriw  fiw 


f^kx{(p  f\(p}l  ^ 


(12) 


and  the  oscillator  strength  due  to  the  z  -  direction  polarization  transition  is 


J  if  • 


ho) 


(•■/Ilk-) 


•  II  \*’/l'S‘k(> 

w  mb\nb^ 


1 

b  2 


^  *  I — 

l__  mw  nw 
V  V  ; 


(13) 

where  n^,  and  n}j  are  the  refractive  index  of  the  well  and  barrier,  respectively. 

It  can  be  seen  from  Eq.(12)  that  in  the  case  of  symmetrical  quantum  well  the  in-plane 
transition  oscillator  strength  will  be  zero.  Denoting  the  ratio  between  the  oscillator 
strength  of  x  polarization  1— >3  transition  and  that  of  z  polarization  1— >2  transition  as  R. 
After  some  manipulations,  we  get 


X\E2-Ed 
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The  results  for  step  MQW  structure  can  be  obtained  in  a  similar  way.  The  results  are 
rather  tedious  and  will  not  be  listed  here. 

The  systems  under  study  are  24Alno35Gao65As/5A  GaAs/SOOAAlosgGaogjAspi], 
55AGaAs/500AAlo.27Gao„As[21],44Alno.34Gao,66As/250AGaAs[15],40Alno.32Gao.68As/500 
AGaAs[34],  40A  Ino3Gao7As/300A  GaAs[31].  The  temperature  used  in  the  calculation  is 
77  K  and  the  parameters  used  in  the  calculation  can  be  found  from  [32].  In  the 
InGaAs/AlGaAs  and  InGaAs/GaAs  system,  lattice  mismatch  between  these  two  systems 
will  cause  strain.  Strain  effects  induced  by  lattice-mismatch  will  modify  the  energy  band 
structure.  For  conduction  band  at  k=0,  strain  will  cause  a  rigid  shift  AEq  .  For  valence 
band  at  k=0,  strain  has  three  effects  on  it:  strain  will  introduce  a  rigid  shift  for  all  the 
three  valence  bands  (heavy-hole,  light-hole  and  split-off  bands),  this  is  caused  by  the 
hydrostatic  component  AEjf,  strain  will  also  cause  the  reduction  of  symmetry  and 
therefore  cause  the  splitting  of  heavy-hole  and  light-hole  bands,  which  is  caused  by  the 
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shear  uniaxial  component  AE^;  strain  can  also  cause  the  coupling  of  the  light-hole  band 
and  the  split-off  band.  AEq  ^//and  zl^^^can  be  expressed,  respectively  as 

A£c-  =  2ac^'‘~-''gi  (16a) 

Cii 

=  (16b) 

Cm 

(16c) 

Cii 

where  Cq  is  the  hydrostatic  deformation  potential  of  conduction  band,  and  by,  are  the 
hydrostatic  and  uniaxial  deformation  potentials  of  valence  bands,  respectively.  Cy ’s  are 
the  elastic  constants,  and  |  is  the  in-plane  strain  component. 

The  energy  band  gaps  under  strain  for  the  heavy  hole,  light  hole  and  split-off 
bands  at  k=0  are  given  respectively  by 

eV"  =  E]  +  ^Ec  +  ^Eu  +  ^Eu  (l^a) 

“  =  £;  +  Afr + ASr  -  A.,  +  \/a5  +  2A»A£,;  +  9AE:')  ( 1 7b) 

=  £>  A£r + A£ A£,;  -  A«  -  ^AI  +  2A„AE,+9aeI]  ( 1  7c) 

where  Aq  and  Eg^  are  the  spin-orbit  splitting  and  energy  band  gap  without  strain, 
respectively. 

As  mentioned  before,  strain  can  cause  the  reduction  of  symmetry.  As  a  result,  the 
effective  mass  of  the  electron  in  the  conduction  band  will  become  anisotropic  (i.e.,  the 
effective  mass  along  the  growth  direction  will  be  different  from  the  in-plane  one).  Under 
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strain,  the  in-  plane  effective  mass  mc\ i  and  the  growth  direction  effective  mass  can 
be  expressed  as  [33] 


Hk: 

ml 


(18a) 


(18b) 


where  mQ  is  the  free  electron  mass,  Ep  is  the  energy  related  to  the  momentum  matrix 
element  of  transition  from  the  valence  band  to  the  lowest  conduction  band,  a  and  /?  are 
quantities  introduced  by  the  strain.  The  variation  of  electron  effective  mass  of  InGaAs  on 
(001)  GaAs  with  In  composition  is  shown  in  figure  1.3.3.  In  our  calculation  we  have 

taken  into  account  this  effect.  We  calculated  the  value  of  R  at  kf=(2mt^\*Ep/^f^^  and  the 
results  are  listed  in  Table  1.1. 


System 

R 

24A  Ino  jsGao  65AS/5  A  GaAs/300A  AIq  3gGao,62As 

7.1x10-' 

55A  GaAs/500A  AIq  j^Ga^yjAs 

3.5x10-' 

44A  Ino34Gao66As/250A  GaAs 

l.lxlO-' 

40A  InojyGao.egAs/SOOA  GaAs 

0.79x10-' 

40  A  Ino  jGao  7AS/3OOA  GaAs 

0.92x10-' 

Table  1.1  Calculated  results  of  /?/  and  R2  for  the  MWIR  and  LWIR  QWIPs  at 
kF=(2m\\*EF/^f^. 
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It  is  noted  from  this  table  that  the  experimentally  observed  normal  incidence 
absorption  cannot  be  explained  by  the  model  depicted  above. 

1 .3.2  The  role  of  strain  on  normal  incidence  absorption 

Several  /t-type  QWIP  structures  showing  strong  normal  incidence  absorption  are 
based  on  strained-layer  InGaAs/AlGaAs  systems  which  have  lattice  mismatch  between 
the  InGaAs  quantum  well  and  GaAs  substrate  and  hence  cause  significant  strain.  We  have 
studied  the  high-strain  InGaAs/AlGaAs  QWIP  to  see  if  the  compressive  strain  is  uniquely 
responsible  for  the  observed  normal  incidence  absorption. 

The  electron-photon  interaction  Hamiltonian  can  be  written  as 

ff=±AWp  (19) 

2c 

where  e  is  the  electronic  charge,  c  is  the  velocity  of  light,  A  is  the  vector  potential,  W  is 
the  inverse  mass  tensor  and  p  is  the  momentum  operator. 

Considering  the  transition  between  the  two  states,  with  the  wave  fiuictions  for  the 
initial  state  denoted  as  y/i  and  the  final  state  denoted  as  y/f,  i.e.. 


(20a) 

(20b) 

It  can  be  shown  that  the  transition  matrix  element  between  the  initial  state  i  and 
final  state/  Mji  is  proportional  to 

A.W„p,+  AyW„p,+  A,A^P,  (21) 

In  the  case  of  biaxial  strain,  the  inverse  mass  tensor  becomes 
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nic 


rric 


mJ 


(22) 


For  the  TE  polarization,  for  example,  taking  A  in  the  x-direction,  since  ffxz-0, 
Mfi  will  be  zero.  This  is  due  to  the  fact  that  in  InGaAs/GaAs  system  the  strain  direction 
coincides  with  the  main  axes  of  energy  ellipsoid,  no  non-diagonal  component  of  inverse 
mass  tensor  will  occur.  It  is  shown  here  that  the  strain  can  not  be  solely  responsible  for 
the  observed  normal  incidence  absorption. 


1.4  Other  possible  normal  incidence  absorption  mechanisms 

1 .4. 1  Interface  roughness 

If  the  interface  between  the  well  and  barrier  is  very  rough,  a  normal  incident  light 
can  be  scattered  into  all  directions  and  cause  the  absorption  of  normal  incident  light. 
However,  modem  MBE  technology  can  grow  very  smooth  interface  in  the  III-V  material 
systems  such  as  GaAs/AlGaAs  and  InGaAs/AlGaAs.  Even  in  the  strained  QWIPs,  since 
the  thickness  of  the  well  is  within  the  critical  thickness  and  the  dislocation  is  localized 
and  hence  the  interface  is  expected  to  be  very  good.  Fluctuations  of  the  order  of  one 
monolayer  always  are  possilbe,  but  it  is  unlikely  that  this  kind  of  fluctuations  can  cause 
normal  incidence  absorption.  Since  the  refractive  indexes  and  effective  masses  of  two 
adjacent  layers  are  usually  quite  similar. 

1 .4.2  The  effect  of  donor  states 
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Recently  Ryzhii  at  al.  proposed  that  doping  might  be  responsible  for  normal 
incidence  absorption  [34].  They  attributed  the  TM  absorption  to  the  transition  from 
quantum  well  ground  state  to  continuum  state  and  TE  absorption  to  the  transition  from 
the  donor  state  to  quantum  well  ground  state  (total  energy  exceeds  the  barrier  height  due 
to  the  in-plane  energy  component).  However,  it  appears  that  their  model  cannot  explain 
the  coincidence  of  the  transition  peaks  of  TE  and  TM  absorption.  Their  model  shows  that 
normal  incidence  does  not  depend  on  the  doping  concentration,  which  is  in  contrast  with 
the  experimental  results  observed  by  Wang  et  al.  [31,35]. 

1.5  Conclusions 

We  have  discussed  theoretically  the  possible  effects  that  might  be  responsible  for 
the  normal  incidence  absorption,  including  spatial  variation  of  the  effective  mass,  the  role 
of  high  order  bands,  strain  effect  etc.  None  of  these  effects  is  able  to  cause  normal 
incidence  absorption  which  is  comparable  to  the  TM  absorption.  Our  calculation  on  the 
QWIPs,  which  exhibit  experimentally  large  normal  incidence  absorption,  shows  that  the 
TE  (normal  incidence)  absorption  is  about  1%  of  the  TM  absorption.  No  intrinsic 
physical  reason  is  responsible  for  the  observed  normal  incidence  absorption.  Therefore, 
the  observed  normal  incidence  absorption  must  be  due  to  some  unknown  experimental 
conditions. 
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Figure  1.2.1(a).  Under  one-  band  constant  electron  effective  mass 
approximation,  in-plane  polarization  (normal 
incidence)  absorption  is  prohibited. 


Figure  1.2.1(b)  Under  one-  band  space- variant  electron  effective  mass 
approximation,  in-plane  polarization  (normal  incidence) 
absorption  is  allowed. 
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Figure  1 .2.2  In  Kane  model,  band  coupling  and  spin-orbit  interaction  cause  the 
spin  degenerate  states  to  form  two  new  degenerate  states  which 
include  both  the  spin-  up  and  spin-  down  states.  Normal  incidence 
absorption  can  occur  between  states  having  same  dominant  spins 
(solid  lines)  and  states  having  different  dominant  spins  (dashed  lines). 
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Figure  1 .2.3  Coupling  of  the  higher  conduction  band  Fu'  with  the  valence  band 
r,5''  (the  coupling  matrix  is  indicated  by  Q'"'')  and  the  lowest 
conduction  band  F,®  (the  coupling  matrix  is  indicated  by  P')  can 
cause  normal  incidence  absorption. 


E, 


E, 


E,  E,  E 


Figure  1 .3.1  Schematic  conduction  band  diagram  of  the  symmetric  and  step 
MQW  structure  under  consideration.  E,  is  the  energy  level  of  the 
bound  state,  E2  and  Ej  are  the  center  energy  levels  of  the  first  and 
second  continuum  states,  respectively. 
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Figure  1.3.2  Coordinates  adopted  in  the  calculation  for  symmetric  (a)  and  step  (b) 
MQW  structure.  In  (a),  is  the  well  width,  is  the  barrier  width  and 
Vb  is  the  barrier  height.  In  (b),  Lw  is  the  well  width,  L^i  and  L^,2  are  the 
first  and  second  barrier  widths,  respectively,  V|,|  and  Vb2  are  the  first 
and  second  barrier  heights,  respectively. 


Figure  1.3.3  Variation  of  the  effective  mass  of  In^jCaj.^As  on  (001)  GaAs. 
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2.  Intersubband  Absorption  Coefficient  in  QWIPs 

The  responsivity  of  a  QWIP  can  be  written  as 

(23) 

where  hvis  the  incident  photon  energy,  rja  is  the  quantum  efficiency,  is  the  escape 
probability,  and  g  is  the  optical  gain.  The  quantum  efficiency  rja  is  given  by 

(24) 

where  I  is  the  length  of  the  high-field  domain,  and  a  is  the  absorption  coefficient.  It  is 
seen  from  Eqs.  (23)  and  (24)  that  absorption  coefficient  plays  an  important  role  in 
determining  the  responsivity  of  QWIPs.  We  have  performed  the  calculation  of  absorption 
coefficients  for  several  QWIP  structures. 

Considering  an  optical  transition  from  the  initial  state  1  to  the  final  state  2,  the 


absorption  coefficient  can  be  written  as 


{E2-E\-fiO})"  +  rn 


where  n  is  the  refractive  index,  co  is  the  incident  light  frequency,  Ej  and  E2  are  the 
energies  of  the  initial  and  final  states,  respectively.  f(E)  is  the  Fermi-Dirac  distribution 
functions.  ri2  is  a  broadening  parameter.  The  momentum  matrix  element  is  given  by 

P6) 

ipl  and  ^2  growth  direction  envelope  functions  for  the  initial  and  final  states, 


respectively. 
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We  have  used  the  above  formula  to  calculate  the  absorption  coefficient  of  some 
QWIPs,  including  two-stack,  two-color  MWIR  QWIP  [21]  and  TC  QWIPs  (  TC  206  and 
TC  301).  In  the  calculation  the  broadening  parameter  was  taken  as  2  meV.  The 
parameters  used  for  calculating  the  absorption  coefficient  of  two-stack  ,  two-color  MWIR 
QWIP  is  listed  in  Table  2.1.  The  parameters  used  for  TC  206  and  TC  310  are  listed  in 
Table  2.2(a)  and  Table  2.2(b),  respectively. 


Table  2. 1 .  Parameters  used  in  the  absorption  coefficient  calculation  of 
a  two-stack,  two-color  MWIR  QWIP  at  300  K. 


GaAs 

m/  (mo) 

0.063 

0.0961 

0.0647 

m,,"  (mo) 

0.063 

0.0961 

0.0578 

Eg  (eV) 

1.424 

1.8979 

1.0629 

n 

3.2989 

3.1308 

3.3693 

Table  2.2(a)  Parameters  used  in  the  absorption  coefficient  calculation 

ofTC206at77  K. 


m/  (mo) 

0.0752 

0.0667 

0.0663 

m,"  (mo) 

0.0752 

0.0622 

0.0643 

Eg  (eV) 

1.6327 

1.238 

1.3717 

n 

3.221 

3.3245 

3.2941 
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Table  2.2(b)  Parameters  used  in  the  absorption  coefficient  calculation  of  TC  301. 
The  temperature  was  taken  as  77  K. 


GaAs 

IiIq  o5G^,95-^S 

0.066 

0.0836 

0.066 

mj'  (mo) 

0.066 

0.0836 

0.0652 

Eg(eV) 

1.508 

1.7574 

1.443 

n 

3.266 

3.18 

3.278 

The  calculated  results  are  shown  in  figure  2.1.  It  is  seen  from  figure  2.1  that  the 
high-strain  QWIP  has  higher  absorption  coefficient  compared  to  that  of  low-strain  QWIP. 
This  is  due  to  several  reasons:  (1)  the  high-strain  QWIP  has  a  deeper  well  and  hence  a 
large  overlap  of  wave  functions  between  the  ground  and  the  excited  states;  (2)  the  doping 
density  in  high-  strain  TC  QWIP  is  higher  and  hence  the  Fermi  energy  is  higher,  which 
can  increase  the  term  of  f(Ej)-f(E2)  in  Eq.  (25)  and  therefore  increase  the  absorption 
coefficient  correspondingly. 

The  calculation  of  absorption  coefficient  shown  here  provides  a  good  basis  for 
further  study  of  the  relationship  among  responsivity,  optical  gain,  and  optical  absorption 
coefficient. 
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400 


(b)  LS  and  HS-TC  QWIP(LWIR). 

Figure  2.1  Absorption  coefficient  of  (a)  high  strain  MWIR  QWIP  and 
(b)  TC  QWIPs. 
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3.  Study  of  Normal  Incidence  Absorption  in  N-type  QWIPs 


3.1  Influence  of  experimental  conditions  on  normal  incidence  absorption 

In  this  section,  we  describe  the  experimental  study  of  the  normal  incidence 
absorption  in  n-type  direct  gap  QWIPs.  All  three  QWIP  samples  discussed  here  were 
grown  on  semi-insulating  GaAs  substrate  (100)  by  molecular  beam  epitaxy  (MBE) 
technique.  Sample- 1  [20]  was  designed  to  have  two  stacks  of  quantum  wells  as  the  active 
region  with  a  contact  layer  in  between  the  two  stacks.  The  bottom  MWIR  stack  has  20 
periods  of  30  nm  Alo.38Gao,62As  barrier  and  2.4  nm  Ino.35Gao.65As  sandwiched  between 
two  0.5  nm  GaAs.  The  top  LWIR  stack  consists  of  20  periods  of  50  nm  Alo.27Gao.73As 
barrier  and  5.5  nm  GaAs  well.  The  Si  doping  density  in  the  wells  is  2.5x10  cm'  for  the 

first  stack  and  0.7x1 0'^cm'^  for  the  second  stack.  Sample-2  has  similar  structure 
parameters  except  that  there  is  only  three  periods  in  each  stack  and  no  contact  layer  in 
between  the  two  stacks.  Sample-3  consists  of  only  three  quantum  wells  with  three 
different  barrier  heights  corresponding  to  three  different  peak  wavelengths  [36].  The 
detector  consists  of  three  5.0  nm  doped  GaAs  wells  and  four  undoped  AlGaAs  barriers. 
There  is  no  InGaAs  involved  and  hence  is  free  from  strain  caused  by  the  lattice  mismatch. 
The  energy  band  diagrams  for  these  three  structures  are  shown  in  figure  3.1.1.  The 
samples  are  processed  using  standard  wet  chemical  etching,  which  produces  certain  edge 
area  on  the  mesa.  The  area  of  the  edge  increases  with  the  etching  depth  of  the  mesa, 
which  depends  on  the  QWIP  structure  thickness.  Sample-1  was  processed  in  two  different 
ways  in  which  the  first  method  was  a  three-step  photolithography.  The  top  stack  was 
etched  down  to  the  middle  contact  which  gives  the  mesa  for  the  LWIR  stack  defined  by 
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the  first  photolithography.  The  MWIR  stack  is  defined  by  etching  the  bottom  stack  to  the 
bottom  contact.  The  third  lithography  defines  the  metalization  area,  which  covers  the 
whole  mesa  area  in  this  processing.  The  second  method  of  processing  was  to  etch  the  first 
stack  down  to  the  middle  contact  for  the  LWIR  stack,  while  the  MWIR  mesa  was  formed 
by  first  totally  removing  the  LWIR  stack  and  then  etching  down  to  the  bottom  contact. 
The  metal  contact  only  occupies  a  small  portion  of  the  mesa  area,  and  front  illumination 
is  used  for  the  photoresponse  measurements.  Sample-2  used  the  second  method  in 
processing.  Sample-3  has  single  color  (LWIR),  the  metalization  used  the  first  processing 
geometry.  The  mesa  areas  are  210x210  pm  for  sample- 1  and  -2  using  second  processing 
method,  and  220x360  pm^  for  sample  -3.  The  ratios  of  the  edge  area  to  the  mesa  area  are 
6.74%,  3.51%  and  0.79%,  respectively,  for  sample  -1,  -2,  and  -3. 

The  photocurrent  spectra  of  the  test  devices  were  measured  using  a 
monochrometer  and  a  calibrated  blackbody  source.  Figure  3.1.2  shows  the  responsivity  of 
sample- 1  as  a  function  of  bias  voltage  under  normal  incidence  without  grating,  along  with 

under  45°  facet  incidence.  The  peak  wavelength  of  the  MWIR  stack  under  normal 
incidence  is  at  4.4  pm  with  a  maximum  responsivity  of  0.2A/W  at  3.2V.  It  can  be  seen 
from  this  figure  that  the  responsivity  ratio  of  the  normal  incidence  to  the  45°  facet 
incidence  is  about  41%  at  2V.  For  the  LWIR  stack,  the  peak  wavelength  of  the  normal 
incidence  is  at  9.2  pm  with  a  maximum  responsivity  of  0.2  A/W  at  2V.  The  responsivity 

ratio  of  the  normal  incidence  to  the  45°  facet  incidence  is  about  40%  at  this  bias.  The  data 
shown  at  2V  used  samples  with  second  processing  method.  Using  first  processing 
method  gave  similar  results.  Figure  3.1.3  shows  the  responsivity  of  sample-2  as  a 
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function  of  bias  voltage  under  normal  incidence  without  grating  and  the  45°  coupling. 
The  peak  wavelength  of  the  MWIR  stack  is  at  4.3  pm  with  a  maximum  responsivity  of 
0.17  A/W  at  2.2  V,  while  the  peak  wavelength  of  the  LWIR  stack  is  at  8.4  pm  with  a 
maximum  responsivity  of  0.1  A/W  at  -3V.  The  responsivity  ratios  of  the  normal 

incidence  to  the  45°  incidence  are  25%  at  2.5  V  and  46%  at  2  V,  for  the  LWIR  and 
MWIR  stack,  respectively.  Figure  3.1.4  shows  the  responsivity  spectra  of  sample-3  for 
the  normal  incidence  without  grating  under  different  biases.  The  responsivity  is  0.49  A/W 
at  the  peak  wavelength  of  7.5  pm  and  bias  voltage  of  -0.58V.  Samples-2  and  -3  have 
much  smaller  edge  areas  compared  with  sample- 1  due  to  their  much  thinner  active  layer 
thickness,  but  they  still  exhibited  a  large  normal  incidence  responsivity. 

The  three  samples  given  above  have  very  different  mesa  heights,  which  are  3358 
nm,  1747  nm  and  505  nm  for  sample- 1,  -2  and  -3,  respectively.  Since  the  top  contact 
layer  also  contributes  to  the  edge  effect,  sample-3  has  not  only  a  very  thin  active  layer, 
but  also  a  100  nm  top  contact  layer  which  gives  a  much  thinner  mesa  height  than  the 
other  two  samples.  The  responsivity  value  of  sample-3  in  Fig.  3.1.4  demonstrates  that  a 
sample  with  a  very  small  edge  area  can  still  have  very  large  normal  incident  responsivity. 
Another  concern  is  the  edge  of  the  substrate,  which  can  also  contribute,  to  the  edge  effect. 
Special  attention  has  been  given  to  eliminate  the  effect  of  the  substrate  edge  by  using  a 
very  small  size  beam  as  shown  in  Fig.  3.1.5.  Normal  incidence  responsivity  was  observed 
in  sample-3  using  an  IR  beam  spot  size  of  317  pm  diameter.  Standard  setup  for  the 
photocurrent  measurements  has  a  blackbody  beam  size  of  around  one  inch  diameter. 
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The  roughness  of  the  metal-semiconductor  interface  at  the  metal  contact  due  to 
the  metal  alloying  inter-difflision  is  another  possibility  that  causes  photon  diffraction  to 
the  in-plane  direction.  Two  metal  contact  geometries  used  to  process  sample- 1  and  -2 
give  two  quite  different  metal  contact  areas.  The  high  responsivity  of  both  sample  -1  and 
-2  reveals  that  metal-semiconductor  interface  diffusion  is  not  the  dominant  factor  for 
normal  incidence  intersubband  transitions. 

Sample-1  and  -2  have  high  concentration  of  indium  in  the  InGaAs/AlGaAs  which 
causes  high  strain  when  grown  on  GaAs  substrate  due  to  the  lattice  mismatch.  The 
material  of  sample-3  uses  only  GaAs/AlGaAs  with  no  InGaAs  involved.  The  high 
responsivity  observed  in  sample-3  at  normal  incidence  demonstrates  that  the  unstrained 
quantum  wells  can  have  high  normal  incidence  absorption  as  well  when  short-periods 
quantum  well  is  used. 

3.2  Investigation  of  new  S-QWIP  structures 

We  designed  and  processed  several  simplified  (S-)  QWIP  structures  for  LWIR 
(QWIP-A)  and  MWIR  (QWIP-B  and  -C)  detection.  We  expect  these  QWIPs  to  show 
large  responsivities  due  to  the  short  period  and  hence  a  relatively  large  photoconductive 
gain.  For  the  MWIR  QWIPs,  a  very  high  indium  composition  of  37%  is  used  in  the 
InGaAs  quantum  well.  The  MWIR  QWIP-C  is  almost  the  same  as  the  MWIR  QWIP-B 
except  that  a  very  small  amount  of  phosphorous  was  added  in  the  AlGaAs  barrier  layer. 
The  phosphorus  is  expected  to  improve  the  interface  property  and  thus  enhances  the 
overall  performance  of  the  device. 
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3.2.1  LWIR  OWIP-A 


The  growth  structure  of  LWIR  QWIP-A  is  shown  in  figure  3.2.1.  The  dark  current 
density  of  this  device  is  shown  in  figure  3.2.2.  The  responsivity  of  this  device  under 
normal  incidence  at  35  K  is  shown  in  figure  3.2.3.  The  responsivity  of  this  device  under 
normal  incidence  at  77  K  is  shown  in  figure  3.2.4.  Figure  3.2.5  shows  the  responsivity  of 
this  device  under  45“  incidence  at  35  K  and  figure  3.2.6  shows  the  responsivity  of  this 
device  under  45“  incidence  at  77  K.  We  have  calculated  the  detectivity  (D*)  value  of  this 
device  under  normal  incidence  at  77  K.  The  results  are  listed  in  Table  3.2.1. 


Table  3.2.1  Detectivity  of  LWIR  QWIP-A  under  normal  incidence  at  77  K. 


bias  (V) 

-1.5 

-1.0 

-0.5 

0.5 

1.0 

1.5 

D‘(lO’cmHz''VW) 

4.17 

3.81 

1.13 

7.33 

3.77 

4.31 

3.2.2  MWIR  OWIP-  B  and  C 

We  have  also  performed  dark  I-V  and  photoresponse  measurements  on  two 
MWIR  QWIPs,  i.e.,  QWIP-B  and  QWIP-C.  The  layer  structure  of  MWIR  QWIP-B  is 
shown  in  figure  3.2.7.  The  layer  structure  of  QWIP-C  is  similar  to  that  of  QWIP-B, 
except  the  barrier  thickness  is  300  A  and  as  mentioned  before,  there  is  a  very  small 
amount  of  phosphorus  added  to  the  barrier  layer.  The  dark  current  density  of  QWIP  467 
and  QWIP  475  are  shown  in  figure  3.2.8  and  figure  3.2.9,  respectively.  For  QWIP-B,  no 
photoresponse  was  observed,  and  the  room  temperature  resistance  of  this  QWIP  device  is 
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always  very  high  (>  1  OK  Q).  For  QWIP-C,  we  were  able  to  measure  responsivity,  and  the 
result  is  shown  in  figure  3.2.10.  The  peak  wavelength  for  this  device  is  at  3.9  pm.  We 
have  calculated  the  detectivity  of  this  device  under  normal  incidence  at  77  K,  and  the 
results  are  listed  in  Table  3.2.2. 


Table  3.2.2  detectivity  of  QWIP  LWIR-C  under  normal  incidence  at  77  K. 


bias  (V) 

-2 

-1.5 

-1 

o 

1 

0.5 

D*  (10*cmHz"'/W) 

0.674 

1.11 

1.55 

1.28 

0.01 

3.3  Normal  incidence  absorption  in  regular  QWIPs  and  EQWIP 

To  explore  the  effect  of  substrate  thinning  on  the  normal  incidence  absorption,  we 
have  performed  a  comparison  study  on  the  following  QWIP  structures: 
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Table  3.3.1  The  layer  structures  for  the  QWIP  under  compared  study. 


Layer 

Thickness  (A) 

Doping  (cm’^) 

GaAs 

4000 

3x10'* 

500 

" 

GaAs 

55 

5x10” 

-^10.27^^.73^^ 

500 

GaAs 

10000 

3x10'* 

1000 

” 

GaAs  Substrate 

- 

The  measurement  of  the  QWIP  device  without  removing  the  substrate  was  done  at 
University  of  Florida  (UF)  and  the  measurement  with  substrate  being  removed  was  done 
at  Lockheed  Martin  Vought  Systems  (LMVS).  The  dark  current  density  of  this  device 
measured  at  UF  is  shown  in  figure  3.3.1.  A  comparison  of  the  dark  current  density 
measured  at  UF  and  LMVS  is  shown  in  table  3.3.2 


Table  3.3.2  Comparison  of  the  dark  current  density  of  QWIP-  D  at  77K. 


Bias  (V) 

Dark  current  density  J  (A/cm^) 

device  V3  of  LMVS 

device  processed  at  UF 

Juf/JlM 

-0.25 

6.78x10' 

1.06x10-' 

1.56 

-0.5 

4.9x10-' 

8.05x10-' 

1.64 

-1.0 

1.9x10-* 

4.3x10-* 

2.26 
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The  peak  responsivity  of  QWIP-D  with  substrate  was  measured  at  both  the  normal 
incidence  and  45°  incidence  illumination,  and  the  results  are  shown  in  figure  3.3.2(a)  and 
figure  3.3.2(b),  respectively.  It  is  seen  from  figure  3.3.2  that  this  device  shows  a  rather 
strong  normal  incidence  absorption.  Table  3.3.3  lists  the  peak  responsivity  of  this  device 
under  normal  incidence  and  45°  incidence. 


Table  3.3.3  Peak  responsivity  of  QWIP-D  under  normal  incidence  and  45° 
incidence  and  their  ratios  at  different  biases. 


Bias  (V) 

normal  incidence 

45°  incidence 

Rp(90°)/Rp(45°) 

Xp(pm) 

Rp(A/W) 

>.p(pm) 

Rp(A/W) 

-0.25 

8.6 

0.039 

8.6 

0.  095 

0.4105 

-0.5 

8.6 

0.209 

8.6 

0.512 

0.4082 

-1.0 

8.6 

0.424 

8.6 

1.05 

0.4038 

The  responsivity  of  this  device  with  substrate  being  removed  was  also  measured 
at  LMVS,  and  no  normal  incidence  absorption  was  found  under  this  condition. 

3.4  Conclusions 

We  have  studied  the  effect  of  experimental  conditions  on  normal  incidence 
absorption  in  an  «-type  direct  gap  QWIP,  including  mesa  form,  mesa  height,  roughness  of 
the  metal-semiconductor  interface,  strain  effect,  substrate  effect  etc.  From  our  study,  we 
have  found  that  the  mesa  form  and  height,  metal-semiconductor  interface  and  strain 
effects  are  not  the  key  factors  which  cause  normal  incidence  absorption  in  this  QWIP. 
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The  samples  used  in  our  study  are  all  grown  on  the  semi-insulating  GaAs  substrate  which 
does  not  absorb  IR  light,  but  scattering  IR  light  diffracted  back  into  the  active  quantum 
well  layers  can  cause  the  observed  normal  incidence  absorption.  The  presence  of 
substrate  could  be  the  key  factor  responsible  for  the  occurrence  of  normal  incidence 
absorption  in  the  n-type  direct  gap  QWIPs. 
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Figure  3.1.1  Schematic  energy  band  diagrams  for  QWIP  samples  -1  (a),  -2  (b)  and  -3(c). 
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RESPONSIVITY(AAV)  RESPONSIVITY(A\W) 


(a)  BIAS  VOLTAGE  (V)  (b)  BIAS  VOLTAGE  (V) 


Figure  3.1.2  Responsivity  of  sample-1  (a)  MWIR  QWIP  and  (b)  LWIR  QWIP. 


(a)  BIAS  VOLTAGE  (V)  (b)  BIAS  VOLTAGE  (V) 


Figure  3.1.3  Responsivity  of  sample-2  (a)  MWIR  QWIP  and  (b)  LWIR  QWIP. 
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Figure  3.1.4  Responsivity  of  sample-3. 
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Figure  3.1.5  Schematic  diagrams  showing  normal  incidence  illumination  for 

(a)  sample-3,  using  backside  illumination  with  a  beam  size  of  3 17pm 
diameter  and  (b)  samples  -1  and  -2,  using  frontside  illumination 
with  a  beam  size  of  1  inch  diameter. 
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Figure  3.2.3  Responsivity  of  LWIR  QWIP-A  under  normal  incident 
illumination  at  35  K. 
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Figure  3.2.4  Responsivity  versus  wavelength  for  the  LWIR  QWIP-A 
with  normal  incidence  illumination  at  77  K. 
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Figure  3.2.5  Responsivity  versus  wavelength  for  the  LWIR  QWIP-A 
with  45°  facet  incidence  illumination  at  35  K. 
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Figure  3.2.6  Responsivity  versus  wavelength  for  the  LWIR  QWIP-A 
with  45°  facet  incidence  illumination  at  77  K. 
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Figure  3.2.8  Dark  current  density  versus  bias  for  the  MWIR 
QWIP-B  measured  at  77  K. 


Figure  3.2.9  Dark  current  density  of  MWIR  QWIP-C  at  77  K. 


Figure  3.2.10  Responsivity  of  MWIR  QWIP-C  under  normal  incidence  at  35  K 
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Figure  3.3.1  Dark  current  density  versus  bias  for  QWIP-  D 
measured  at  77  K  (with  no  substrate  thinning). 


Figure  3.3.2  (a)  Responsivity  versus  wavelength  for  QWIP-  D  with 

normal  incidence  illumination  measured  at  different  biases. 
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Figure  3.3.2  (b)  Responsivity  versus  wavelength  for  QWIP-  D  with 
45°  incidence  illumination  at  different  biases. 
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4.  Simplified  (S-)  QWIPs  for  Low  Temperature  and  Low  Background 

Applications 


4.1  Variation  of  Black  Body  Flux  With  Field  of  View  Angle 

One  important  application  area  for  QWIPs  is  in  the  situation  of  low  temperature 

and  low  background  environments.  It  is  important  to  study  the  variation  of  blackbody 
flux  with  field  of  view  (FOV)  angle.  The  photon  flux  (j)  within  wavelength  range  of  X,,  and 
X2  can  be  formulated  as 


exp(  ^  -  I 


.  2^ 
Sin  — 


(27) 


^ksT 


where  Q  is  the  FOV  angle  (see  figure  4.1).  We  have  calculated  the  photon  flux  within  0- 
100  pm  at  different  temperatures,  and  the  results  are  shown  in  figure  4.2. 


4.2  Low  Temperature  and  Low  Background  Modeling  for  the  S-QWIP 

The  simplified  (S-)  QWIPs  are  QWIPs  with  reduced  number  of  quantum  wells 
(typically  less  than  five)  and  they  offer  some  unique  advantages,  such  as  shorter  growth 
time,  higher  yield,  higher  array  uniformity  etc.  At  low  background  the  device  has  to  be 
operated  at  lower  temperatures,  and  it  is  difficult  to  achieve  BLIP  operation  at  a  low 
background  even  with  reduced  operation  temperature.  It  is  useful  to  predict  theoretically 
at  what  background  the  device  will  achieve  BLIP  operation.  Here  we  consider  sample  3 
described  in  section  3.1.  Figure  4.3  shows  the  dark  current  of  this  device  at  77K  and  40K, 
as  well  as  the  300K  window  current  when  the  detector  is  kept  at  77K.  The  dashed  lines  in 
figure  4.3  gives  the  calculated  window  current  with  a  photon  flux  of  7.7x10”  cm'^s'', 
4.5x10’'*  cm'^s  '  and  l.lxlO’^cm'^s'’,  from  top  to  the  bottom,  respectively.  From  this  figure 
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one  can  see  that  the  device  is  BLIP  at  77K  with  a  300K  background  and  36°  field  of  view 
angle.  The  top  dash  line  is  the  calculated  300K  background  current  which  matches  well 
with  the  measured  data.  When  the  photon  flux  is  reduced,  the  device  is  no  longer  BLIP  at 
77K.  From  the  figure  one  can  see  that  the  device  is  BLIP  at  40K  with  a  low  background 
of  4.5x1 0'‘’cm'^s''.  When  the  FOV  is  reduced  as  in  most  of  the  systems,  the  background 
photon  flux  has  to  be  higher  in  order  for  the  detector  to  have  BLIP  operation. 

At  low  background  with  a  far  away  target,  thermal  imaging  is  very  difficult  and 
we  usually  see  an  unresolved  target  with  an  IR  FPA.  An  infi-ared  system  model  was  used 
to  calculate  the  target  detection  range  as  a  function  of  several  detector  parameters.  Path 
transmission  was  assumed  to  be  1  for  simplicity.  A  dark  current  density  of  10’’  A/cm^ 
was  used  for  a  detector  temperature  of  40  K,  and  10'^  A/cm^  was  used  for  77  K  detector 
temperature.  The  conversion  efficiency  is  defined  as  photon  absorption  quantum 
efficiency  times  photoconductive  gain.  The  quantum  efficiency  is  assumed  to  be  25% 
while  the  gain  varies  to  change  the  conversion  efficiency.  At  77K,  the  photoconductive 
gain  is  assumed  to  be  equal  to  dark  current  noise  gain.  The  integration  time  is  varied  to 
fill  half  of  the  charge  wells  of  the  ROIC.  At  40K,  the  dark  current  noise  gain  is  assumed 
to  be  0.2  which  is  not  equal  to  photoconductive  gain.  The  integration  time  is  set  to  be  30 
ms  and  the  charge  wells  are  not  quite  filled  to  half  of  the  wells.  The  parameters  which  do 
not  change  in  the  modeling  are  given  in  Table  4.2.1. 
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Table  4.2.1  Parameters  used  in  the  modeling 


Parameters 

Value 

Units 

Target  Size 

1x1 

Background  Flux 

Window  Temperature 

40 

K 

Optics  Temperature 

40 

K 

Cold  Shield  Temperature 

40 

K 

Peak  Wavelength 

8.6 

|im 

Spectral  Bandwidth 

1.5 

|im 

Pixel  Size  /  Pixel  Pitch 

38/40 

jam 

Noise  Gain 

0.2 

Charge  Well  Capacity 

o 

X 

Electrons 

False  Alarm  Rate 

IxlO-'^" 

Signal  to  Noise  Ratio 

4.8 

Probability  of  Detection 

50 

% 

The  background  flux  is  chosen  to  be  1x10'^  photons/cmVsec  which  is  not  in  BLIP 
condition.  When  the  target  is  maintained  at  the  same  temperatures,  lower  background 
photon  flux  will  give  a  longer  detection  range.  For  40  K  operation,  the  dark  current  and 
the  multiplexer  (MUX)  noises  are  about  half/half,  assuming  the  MUX  noise  is  300 
electrons.  At  77  K,  the  dark  current  noise  dominates.  The  target  temperatures  are  set  at 
240K,  270K  and  350  K.  The  results  are  given  in  figure  4.4  in  which  the  detection  range  is 
plotted  as  a  function  of  conversion  efficiency.  The  unresolved  target  detection  range  is 
defined  as  having  a  signal  to  noise  ratio  of  4.8  and  probability  of  detection  50%.  A 
regular  QWIP  gives  a  conversion  efficiency  between  2  to  6  %.  When  the  conversion 
efficiency  increases,  the  detection  range  increases  as  seen  in  the  figure  2.  At  40  K, 
assuming  the  dark  current  noise  gain  is  different  from  the  photoconductive  gain,  a  more 
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significant  increase  of  detection  range  is  achieved  than  at  77  K  where  dark  current  gain  is 
similar  to  photoconductive  gain. 

The  optics  has  been  assumed  to  be  cooled  at  40K  which  does  not  contribute  to  the 
total  noise  in  the  modeling.  If  the  sensor  system  is  not  in  space  and  the  optics  is  warm  at 
300K,  40K  operation  will  have  optics  noise  dominate.  The  detection  range  will  decrease 
to  similar  to  77K  operation  due  to  the  optics  noise.  Special  readout  design  to  subtract  the 
warm  optics  current  is  necessary  for  long  range  detection.  Another  factor  is  that  we  have 
assumed  the  path  transmission  as  1  in  all  cases.  When  the  background  is  high,  the  path 
transmission  is  smaller  than  1.  Figure  4.5  gives  a  basic  idea  how  the  path  transmission 
coefficient  can  dramatically  change  the  detection  range. 
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Field  of  View  Angle  6 


Figure  4.1  Schematic  illustration  of  field  of  view  angle. 


FIELD  OF  VIEW  ANGLE  (n) 


Figure  4.2  Variation  of  photon  flux  (photons  per  unit  area  per  unit  time)  with 

field  of  view  angle  at  different  temperatures.  The  wavelength  range 
under  consideration  is  from  0  to  100  pm. 
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Figure  4.3  The  dark  current  and  window  photocurrent  of  a  three-well  S-QWIP 
with  calculated  window  photocurrent  at  different  background 
photon  fluxes. 
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Figure  2.  The  detection  range  vs.  caiversion  efficiency  at  target  tenperatures  of 
240K,270Kand350K 

Figure  4.4 
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Figure  4.5  The  detection  range  as  a  function  of  path  transmission. 
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5,  A  Novel  Four-  Stack  Four-  Color  QWIP  for  MWIR  and  LWIR  Applications 


5.1  Introduction 

Recently  there  is  an  increasing  interest  in  the  development  of  multi-color  QWIPs 
for  MWIR  and  LWIR  dual  band  detection.  Development  of  a  QWIP  focal  plane  array 
(FPA)  for  staring  image  sensor  with  both  MWIR  and  LWIR  spectral  bands  on  single  chip 
is  a  very  important  step  in  meeting  the  needs  and  technological  challenges  of  future 
imaging  applications,  including  space  exploration  and  sensing  systems,  medical  imaging 
and  basic  scientific  research.  Two  common  approaches  have  been  used  to  achieve  multi¬ 
color  detection.  One  common  approach  is  to  employ  coupled  quantum  wells  which  is 
voltage-tunable  to  realize  multi-color  detection  [37,38].  However,  the  wavelength 
tunability  is  usually  quite  small.  Another  approach  uses  the  stacked  structure  consisting 
of  two  (or  more)  stacks  of  QWIPs  with  one  MWIR  stack  and  one  LWIR  stack  grown  on 
the  semi-insulating  GaAs  substrate.  In  this  way  one  can  realize  simultaneous  detection  of 
MWIR  and  LWIR  wavebands.  Employing  this  type  of  structure,  Tidrow  et  al.  have  made 
a  high  quality  two-stack,  two-color  QWIP  for  MWIR  and  LWIR  detection[20]. 

The  simplified  QWIP  (S-QWIP)  discussed  earlier  usually  contains  quantum  wells 
of  less  than  ten  periods,  whereas  the  conventional  QWIPs  contain  40  to  50  periods  of 
active  quantum  wells.  The  S-QWIP  has  many  advantages,  such  as  reduced  material 
growth  time  and  high  yield,  easier  array  fabrication,  low  bias  and  hence  lower  power 
consumption,  easier  to  realize  multi-color  in  a  stacked  structure  etc.  These  advantages 
have  made  the  S-QWIP  very  attractive  for  multi-color  multi-stack  QWIP  applications.  A 
three- well  QWIP  has  been  demonstrated  which  shows  very  high  performance  [36].  A 
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very  unique  feature  about  S-QWIP  is  that  the  optical  gain  of  it  is  very  large  compared  to 
that  of  long  -  period  QWIPs  due  to  the  reduced  transit  time  of  photo-excited  electrons. 
This  makes  S-QWIPs  very  useful  for  low  photon  background  detection. 

The  schematic  energy  band  diagram  of  the  4-color  QWIP  device  is  shown  in 
figure  5.1.  It  consists  of  two  stacks  of  multiple  quantum  wells  as  the  active  region  with  a 
highly  doped  contact  layer  in  between  them.  The  first  stack  is  designed  to  detect  one 
LWIR  and  one  MWIR,  the  LWIR  structure  consists  of  eight  periods  of  a  600  A 
Alo.27Gao,73As  barrier  and  a  55  A  GaAs  well,  the  doping  density  is  0.7x1 0'*cm'^  the 
MWIR  structure  consists  of  four  periods  of  a  400  A  AlojgGao^jAs  barrier  and  a  35  A 
Ino22Gao7gAs  well,  the  doping  density  is  2.8x10'*  cm'*.  The  second  stack  is  designed  to 
detect  two  different  LWIRs,  the  first  LWIR  structure  consists  of  five  periods  of  a  500  A 
Alo.25Gao 75AS  barrier  and  a  50  A  GaAs  well,  the  doping  density  is  0.6x10'*  cm*;  the 
second  LWIR  structure  consists  of  ten  periods  of  a  800  A  AIq  ,6Gao  g4As  barrier  and  a  65 
A  GaAs  well,  the  doping  density  is  0.4x10'*  cm'*.  The  middle  contact  is  made  of  4000  A 
GaAs  and  the  doping  density  is  1x10'*  cm'*.  The  sample  was  grown  on  semi-insulating 
(100)  GaAs  by  molecular  beam  epitaxy  (MBE)  using  Asj  as  the  arsenic  source. 

The  sample  was  processed  into  devices  with  area  of  216x216  pm*  using  standard 
photolithography  and  wet  chemical  etching.  The  mesa  structures  of  the  top  stack  and 
whole  stack  were  formed  by  etching  down  from  the  top  contact  layer  to  the  middle 
contact  layer  and  to  the  bottom  contact  layer,  respectively.  The  mesa  structure  of  the 
bottom  stack  was  formed  first  by  completely  removing  the  top  contact  and  the  top  stack, 
and  then  etching  down  from  the  middle  contact  layer  to  the  bottom  contact  layer.  Ohmic 
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contact  was  made  by  depositing  a  square  contact  ring  composed  of  AuGe/Ni/Au  on  the 
periphery  of  the  mesa.  The  device  substrate  was  polished  to  have  a  45°  facet  for  IR 
illumination.  No  substrate  thinning  was  made  on  this  device.  In  the  dark  current  and 
photoresponse  measurements,  the  bottom  contact  layer  was  grounded  and  top  (middle) 
contact  layers  were  positively  or  negatively  biased.  The  photoresponse  was  measured  by 
back  illumination  using  a  monochrometer  and  a  calibrated  blackbody  source. 

5.2  Device  Performance  and  Characteristics 


5.2.1  The  top  stack 

Figure  5.2  shows  the  dark  current  density  of  the  top  stack.  The  solid  lines  are  dark 
current  densities  measured  at  40K,  60K,  77K  and  90K  and  the  dashed  line  is  the 
measured  300K  window  current  density  with  a  180°  FOV.  The  BLIP  temperature  was 
found  to  be  60K  up  to  ±2.5V. 

Figure  5.3(a)  shows  the  responsivity  of  top  stack  at  40K  under  several  different 
bias.  Two  peaks  were  observed,  one  at  9.0  |j,m  and  another  one  at  12.3  pm.  At  ±1.5V  and 
±1V,  two  peaks  appear  simultaneously  while  at  +3V  only  the  12.3  pm  peak  shows  up 
with  a  responsivity  larger  than  2  A/W.  The  full  width  at  half-maximum  (FWHM)  for  both 
peaks  are  around  1.3  pm.  Figure  5.3(b)  shows  the  peak  responsivities  at  9.0  pm  and  12.3 
pm  versus  bias  voltage  at  40K.  The  responsivities  at  both  peak  wavelengths  increase 
initially  with  applied  voltage  and  after  reaching  a  maximum  value,  the  responsivity  with 
peak  wavelength  at  9.0  pm  drops  continually.  The  responsivity  with  peak  wavelength  at 
12.3  pm  first  reaches  a  maximum  value,  then  drops  continually,  after  reaching  a 
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minimum  value  it  increases  again.  One  can  see  from  this  figure  that  the  two  peak 
responsivities  exhibit  an  interesting  competition  with  applied  voltage. 

When  the  temperature  is  increased  to  77K,  the  12.3  pm  response  peak  is  almost 
not  observable,  as  can  be  seen  from  figure  5.4(a)  which  shows  the  peak  responsivities 
versus  applied  voltage.  The  is  due  to  the  fact  that  at  77K  the  thermionic  emission  dark 
current  for  the  12.3  pm  peak  is  too  large,  which  prevents  the  photocurrent  to  be  observed 
at  this  temperature.  Figure  5.4(b)  shows  the  spectral  responsivity  of  the  top  stack  at  0.5V 
and  -0.5V  bias  and  77  K.  As  can  be  seen  from  this  figure,  only  the  9.0  pm  peak  appears. 
Compared  to  the  40K  spectral  responsivity,  the  peak  position  does  not  change  but  the 
responsivity  increases  at  77K. 


5.2.2  The  bottom  stack 

Figure  5.5  shows  the  dark  current  density  of  the  bottom  stack  at  several  different 
temperatures  along  with  the  300K  window  current  density  with  180°  FOV.  The  BLIP 
temperature  was  found  to  be  77K  up  to  ±5  V. 

Figure  5.6(a)  shows  the  LWIR  responsivity  of  bottom  stack  at  3V  and  -4V  bias, 
and  77K.  The  peak  wavelength  is  found  to  be  at  8.5  pm  under  both  biases.  At  3V,  the  cut- 
on  full  width  at  half-maximum  (FWHM)  is  at  8.1  pm  and  cutoff  FWHM  at  9.3  pm  with  a 
spectral  bandwidth  of  1.2  pm.  At  -4V,  the  cut-on  FWHM  is  at  8.1  pm  and  cutoff  FWHM 
at  9.0  pm  with  a  spectral  bandwidth  of  0.9  pm.  Figure  5.6(b)  shows  the  dependence  of 
peak  responsivity  on  the  applied  voltage.  As  can  be  seen  from  this  figure,  the  peak 
responsivity  increases  with  bias.  We  have  also  measured  the  LWIR  responsivity  of  the 
bottom  stack  at  77K,  they  show  a  similar  result  as  at  40K.  Figure5.7(a)  shows  the  LWIR 
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spectral  responsivity  of  the  bottom  stack  at  77K  with  bias  of  3V  and  -4V.  Figure  5.7(b) 
shows  the  bias  dependence  of  the  peak  responsivity.  Compared  to  the  40K  case,  there  is 
an  increase  in  the  responsivity  at  77  K. 

Figure  5.8(a)  shows  the  MWIR  responsivity  of  the  bottom  stack  at  77K  under  two 
different  biases,  2V  and  -2V.  The  response  peak  for  this  MWIR  stack  is  at  4.7  pm.  The 
FWHM  is  about  0.5  pm.  Figure  5.8(b)  shows  the  bias  dependence  of  peak  responsivity. 
As  can  be  seen  from  this  figure,  the  peak  responsivity  increases  with  bias  up  to  +2V 
under  positive  bias  and  -2.5V  under  negative  bias,  and  becomes  saturated  with  a 
maximum  value  of  0.46  A/W  at  2V.  The  spectral  responsivities  at  ±  2V  bias  at  90K  are 
shown  in  figure  5.9(a),  and  the  peak  responsivity  dependence  on  bias  at  90K  is  shown  in 
figure  5.9(b),  respectively.  At  90K,  the  responsivity  under  negative  bias  is  slightly  larger 
than  that  under  the  corresponding  positive  bias. 

5.2.3  The  MWIR-LWIR  stacked  OWIP 

The  dark  current  density  of  the  whole  MWIR-LWIR  stacked  is  shovm  in  figure 
5.10.  As  can  be  seen  from  this  figure,  the  device  is  under  BLIP  at  77K  with  bias  up  to 
±7V.  Figure  5.11(a)  shows  the  LWIR  responsivity  of  the  whole  stack  at  40K  under 
several  different  biases.  Two  response  peaks  at  8.8  pm  and  12.3  pm  were  observed  in  the 
whole  stack  QWIP.  FWHM  for  both  peaks  are  around  1 .3  pm.  It  is  noted  that  there  are 
two  LWIR  QWIPs  contributed  to  the  8.8  pm  peak,  one  at  the  top  stack  with  peak 
wavelength  at  9.0  pm  and  the  other  at  the  bottom  stack  with  peak  wavelength  at  8.5  pm. 
The  interaction  between  these  two  peaks  at  a  certain  bias  results  in  the  observed  8.8  pm 
peak.  Figure  5.11(b)  shows  the  peak  responsivity  versus  applied  bias.  As  can  be  seen 
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from  this  figure,  they  exhibit  a  similar  competition  phenomena  as  we  have  already  seen  in 
the  top  stack  case.  Also  similar  to  the  top  stack  case,  when  the  temperature  is  increased, 
the  12.3  pm  peak  cannot  be  observed  due  to  the  very  large  thermionic  emission  dark 
current  at  high  temperatures.  Figure  5.12(a)  shows  the  LWIR  responsivity  of  the  top  stack 
at  77K  under  two  different  biases,  7V  and  -6V.  Figure  5.12(b)  shows  the  bias  dependence 
of  the  8.8  pm  peak  responsivity  at  77K. 

Figure  5.13(a)  shows  the  MWIR  responsivity  of  the  whole  stack  at  77K  under  two 
different  biases,  3V  and  -6V.  The  peak  wavelength  is  also  at  4.7  pm.  At  both  biases  the 
cut-on  FWHM  is  at  4.35  pm  and  the  cutoff  FWHM  is  at  4.95  pm  with  a  spectral 
bandwidth  of  0.9  pm.  The  dependence  of  peak  responsivity  on  bias  at  77K  is  shown  in 
figure  5.13(b).  The  peak  responsivity  increases  linearly  with  bias  up  to  ±3  V  and  becomes 
saturated  with  a  maximum  value  of  1  AAV  at  3V,  Figure  5.14(a)  shows  the  MWIR 
spectral  responsivity  of  the  whole  stack  at  2V  and  3V.  Figure  5.14(b)  shows  the  bias 
dependence  of  the  peak  responsivity  at  90K.  They  show  a  similar  characteristics  as  the 
77K  case. 
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Figure  5.1  Schematic  energy  diagram  of  the  four-stack,  four-color  QWIP 
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Figure  5.2  Dark  current  density  versus  bias  for  the  top  stack. 
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Figure  5.3(a)  Responsivity  of  the  top  stack  at  40K. 
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Figure  5.3(b)  Peak  responsivity  versus  applied  voltage  for 
the  top  stack  at  40K 
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Figure  5.4(a)  Peak  responsivity  versus  applied  voltage  for 
the  top  stack  at  77K. 
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Figure  5.4(b)  Responsivity  of  the  top  stack  at  77K. 
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Figure  5.5  Dark  current  density  of  the  bottom  stack. 


Figure  5.6(a)  LWIR  responsivity  of  the  bottom  stack  at  40K. 
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Figure  5.6(b)  LWIR  peak  responsivity  of  the  bottom  stack  at  40K. 
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Figure  5.7(a)  LWIR  responsivity  of  the  bottom  stack  at  77K. 


62 


PEAKRESPONSIVITY(WW) 


-4  -3  .2  .1  0  1  2  3  4 


B  IAS  VO  LTAG  E  (V) 


Figure  5.7(b)  LWIR  peak  responsivity  versus  applied  voltage  at  77K 


Figure  5.8  (a)  MWIR  responsivity  of  the  bottom  stack  at  77K. 
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Figure  5.8(b)  MWIR  peak  responsivity  versus  applied  voltage  for 
the  bottom  stack  at  77K. 
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Figure  5.9(a)  MWIR  responsivity  of  the  bottom  stack  at  90K. 
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5.9(b)  MWIR  peak  responsivity  versus  applied  voltage  at  90K. 
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Figure  5.10  Dark  current  density  of  the  whole  stack. 
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Figure  5.11  (a)  LWIR  responsivity  of  the  whole  stack  at  40K. 
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Figure  5.1 1  (b)  LWIR  peak  responsivity  versus  applied  voltage  for  the 
whole  stack  at  40K. 
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Figure  5.12(a)  LWIR  responsivity  of  the  whole  stack  at  77K. 
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Figure  5.12  (b)  LWIR  peak  responsivity  versus  applied  voltage 
for  the  whole  stack  at  77K. 
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Figure  5.13(a)  MWIR  responsivity  of  whole  stack  at  77K. 
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Figure  5.13(b)  MWIR  peak  responsivity  versus  applied  voltage 
of  the  whole  stack  at  77K. 
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Figure  5.14(a)  MWIR  responsivity  of  the  whole  stack  at  90K. 
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Figure  5.14(b)  MWIR  peak  responsivity  versus  applied  voltage  at  90K. 
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6.  Quantum  Confined  Stark  Effect  in  an  Asymmetric  Weakly  Coupled  Double 

Quantum  Wells  Structure 

6.1  Introduction 

Two  common  approaches  have  been  used  to  achieve  multicolor  detection.  One 
approach  uses  the  stacked  structure  consisting  of  two  (or  multi)  stacks  of  QWIPs  with  one 
MWIR  stack  and  one  LWIR  stack  grown  on  the  SI  GaAs  substrate.  Employing  this  kind 
of  structure,  Tidrow  et  al.  made  a  high  quality  two-stack,  two-color  QWIP  [21].  Another 
common  approach  is  to  employ  the  coupled  QW  structure  which  is  voltage-tunable  to 
realize  multicolor  detection  [37,38].  This  wavelength  tunability  is  very  useful  in  military 
application  such  as  target  discrimination.  From  the  point  of  view  of  physics,  voltage¬ 
tuning  is  associated  with  the  motion  of  confined  electrons  and  holes  under  an  electric 
field  along  the  confined  direction.  In  this  case  the  quantum  confined  Stark  effect  is  a  very 
important  issue.  Study  on  the  quantum  Stark  effect  will  help  to  design  new  novel 
multicolor  QWIPs. 

In  general  Stark  shift  refers  to  the  shift  of  energy  levels  of  a  charge  system  under  an 
applied  electric  field.  There  have  been  extensive  study  on  quantum  confined  Stark  effect 
in  quantum  wells  [39-47].  Stark  shift  in  a  single  quantum  well  is  very  small  [39].  To 
increase  the  Stark  shift,  a  high  electric  field  is  needed  and  this  will  induce  large  leakage 
current  in  optical  modulators  and  dark  current  in  QWIPs.  Stark  shift  in  a  step  quantum 
well  structure  has  also  been  studied  [40]  and  it  shows  a  very  large  Stark  shift  compared  to 
the  single  quantum  well  under  the  same  bias.  Other  structures  have  also  been  suggested  to 
further  enhance  Stark  shift.  It  is  very  desirable  to  have  a  large  Stark  shift  under  low  bias 
for  the  practical  applications  of  the  Stark  shift  QWIP  structures. 
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In  this  section  we  describe  an  asymmetric  weakly  coupled  double  quantum  wells 
structure  which  can  have  very  large  Stark  shift  under  very  low  bias.  As  schematically 
shown  in  figure  6.1,  it  consists  of  a  shallow  well  and  a  deep  well  separated  by  a  thin 
barrier  and  confined  by  two  thick  barriers.  The  barriers  are  consisted  of  AloaGa^jAs  and 
the  shallow  well  is  composed  of  Alo  jGao  gAs.  The  deep  well  is  composed  of  GaAs  and 
doping  is  made  in  this  well  with  a  doping  density  of  7x10”  cm‘^.  In  our  discussion  the 
conduction  band  offset  parameter  is  taken  as  0.55  and  all  the  parameters  used  in  the 
calculation  can  be  found  in  [32]. 


6.2  Formulation 

The  wave  function  for  an  electron  moving  in  a  quantum  well  can  be  written  as 

'!'{•■) ''4^)  (28) 

where  S  is  the  area  of  the  well,  k  is  the  wave  vector  of  the  electron  in  the  x-y  plane,  r|  |  is 
the  position  vector  in  the  x-y  plane  and  Uc(r)  is  the  cell  periodic  function  near  the 
conduction-band  extremum,  (^z)  is  the  envelope  function  in  the  z-  (growth)  direction.  The 
Schrodinger  equation  for  in  the  presence  of  an  electric  field  F  (as  shown  in  Fig. 
6.1(b))  is 

■y  r  \ 

+  {V  -eFz)(l>{z)=E(l>{z)  (29) 

where  mj*  is  the  effective  mass  of  electron  in  the  j-th  layer,  Vj  is  the  potential  in  y-th 
layer,  E  is  the  energy  of  electron  associated  with  the  motion  in  z-  direction.  The  double 
wells  system  can  be  biased  in  two  different  ways.  For  the  electric  field  shown  in  figure 
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6.1(b)  we  call  it  under  negative  bias.  In  each  layer  mj*  is  a  constant  and  has  a  sharp 
change  at  the  interface.  The  variation  of  effective  at  the  interfaces  has  a  significant  effect 
on  the  electronic  properties  of  the  quantum  well  system.  Within  each  layer  Eq.  (29)  can 
be  rewritten  as 


ti  ^(z) 


(Vj  -  eFz)0^2)  =  £^(z) 


2Wy  dz 

By  transforming  into  a  dimensionless  coordinate  Z, 


Z  =  - 
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(30) 


(31) 


Equation  (30)  becomes 

-^^z)-z<p{z)  =  o 
dz 

The  solution  of  this  equation  can  be  written  as 


(32) 


(33) 


where  Aj  and  Bj  are  the  Airy  functions  [48].  The  boundary  conditions  are  that  and 
mj*d  d<f)(z)/dz  are  continuous  at  each  interface. 

When  calculating  the  positions  of  the  energy  levels,  the  exchange  interaction 
energy  in  the  ground  state  should  be  taken  into  account,  which  can  be  expressed  as  [49] 


e  dk  k 
0^  {27ty 


(34) 


where  cr  is  the  two-dimensional  electron  density  in  the  quantum  well,  fis 

the  dielectric  constant,  k  and  k'  are  the  in-  plane  wave  vectors  of  electrons.  Under  usual 
doping  condition,  the  exchange  energy  will  lower  the  ground  state  subband  energy  by  5 
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to  10  meV  and  shift  the  absorption  peak  to  shorter  wavelength.  The  direct  interband 
Coulomb  interaction  (repulsive)  is  much  smaller  than  the  exchange  energy  and  can  be 
neglected. 

Considering  an  optical  transition  from  initial  state  1  to  final  state  2  in  the  quantum 
well,  the  oscillator  strength  /  associated  with  this  transition  is 


(35) 


where  my^/*  is  the  effective  mass  of  electron  in  the  well  and  pz  is  the  momentum  operator 
of  electron  in  the  z  -  direction.  The  absorption  coefficient  of  a  system  is  proportional  to 
the  oscillator  strength  and  joint  density  of  states.  Oscillator  strength  plays  an  important 
role  in  determining  the  optical  properties  of  a  system. 

The  absorption  coefficient  can  be  obtained  from  Fermi’s  golden  rule  and  the  result 
is 


(36) 


where  n  is  the  refractive  index,  oo  is  the  incident  light  frequency,  Ej  and  E2  are  the 
energies  of  the  initial  and  final  states,  respectively.  f(E)  is  the  Fermi-Dirac  distribution 
functions.  772  is  a  broadening  parameter.  The  momentum  matrix  element /?72  is 


We  will  use  the  following  formula  for  the  Fermi  energy 

E'=k»T{^ 
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where  k]g  is  the  Boltzmann’s  constant,  Nj;)  is  the  doping  density  and 
n(f=^%m-w*kBTIlfLyv,  T  is  the  temperature,  h  is  the  Planck’s  constant  and  Ly^  is  the  well 
width. 

We  next  consider  the  device  performance  of  this  structure.  For  practical 
applications  of  quantum  well  system  for  QWIPs,  we  are  concerned  with  the  dark  current 
and  responsivity  of  the  structure.  In  order  to  calculate  the  dark  current,  we  need  to 
consider  the  turmeling  probability  of  a  quantum  well.  Assuming  that  the  well  width  is  Ly; 
and  the  barrier  width  is  L}j,  the  effective  mass  of  electron  is  iriyy*  in  well  and  in 
barrier  and  the  barrier  height  is  E}),  then  under  a  bias  F,  the  tunneling  probability  can  be 
obtained  from  the  WKB  approximation  and  the  results  yields 


t{e,f)  = 


exp  - 


3heF 


{2m^'\Eb-E-eF^ 


exp  - 
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E>  Et,-eF  Z,„ 
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(39) 


Under  negative  bias,  the  dark  current  density  can  be  obtained  as 


J,= 


]f(E)ZdE+^2^ — , 


]fiE)T,T^ciE 

E, 


(40) 


where  and  rriy/i*  are  the  effective  masses  of  electrons  in  the  right  well  and  left  well, 
respectively,  Hr  and  jui  are  the  mobility  for  electrons  in  the  right  and  left  wells, 
respectively,  and  v_y/  are  their  respective  saturated  velocities.  Lp  is  the  thickness  of  the 
whole  structure.  Tf  and  F/  are  the  tunneling  probabilities  for  the  right  and  left  wells, 
respectively. 


74 


The  peak  responsivity  Rp  is  given  by 

(41) 

where  hv  is  the  energy  of  incident  light,  a  is  the  absorption  coefficient  of  the  QWIP,  I  is 
the  length  of  high  field  domain  in  QWIP,  pg  is  the  escape  probability  of  electrons  from 
the  well  and  g  is  the  optical  gain. 

We  will  use  the  formula  developed  in  this  section  to  study  the  Stark  effect  in  the 
structure  shown  in  figure  6.1  and  discussed  its  device  performance. 

6.3  Results  and  Discussion 

In  the  following  discussion,  we  shall  first  show  how  the  structure  parameters 
affect  the  energy  levels  and  oscillator  strength,  then  we  will  choose  a  particular  parameter 
and  discuss  the  Stark  effect  and  finally  we  will  discuss  the  device  performance. 

6.3.1  Variation  of  energy  levels  and  oscillators  with  structure  parameters 

There  are  two  bound  levels  in  the  proposed  QWIP  structure,  the  lower  level  is 
from  the  deep  well,  and  the  higher  level  is  from  the  shallow  well.  The  position  of  these 
two  levels  will  be  affected  by  the  structure  parameters,  such  as  the  left  well  width,  thin 
barrier  width  and  right  well  width.  In  our  discussion  the  thick  barrier  thickness  is  kept  at 
300  A. 

Figure  6.2  shows  the  effect  of  right  well  width  on  the  E ]  eind  E2  and  the  oscillator 
strength  associated  with  the  transition  between  the  two  levels.  In  the  calculation,  the  left 
well  width  is  set  to  be  70  A,  and  the  thin  barrier  thickness  is  taken  to  be  30  A.  It  is  seen 
from  figure  6.2(a)  that  Ei  decreases  rapidly  with  right  well  width,  while  E2  decreases 
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slowly.  This  can  be  understood  from  a  single  quantum  well  case.  It  is  well  known  that  in 
a  single  quantum  well,  the  energy  levels  decrease  with  the  well  width.  Since  Ej  is 
originated  from  the  right  well  which  is  weakly  coupled  to  the  left  well,  it  will  be  greatly 
affected  by  the  right  well  width.  Due  to  the  coupling  with  the  right  well,  the  effective  left 
well  width  will  be  increased  with  the  right  well  width,  and  the  E2  level  which  is 
originated  from  the  left  well  will  also  decrease  correspondingly.  It  is  also  seen  from 
figure  6.2(a)  that  exchange  energy  will  cause  the  ground  state  (Ej)  to  be  lowered  down 
about  1 0  meV.  The  oscillator  strength  first  decreases  with  the  right  well  width  and  then 
almost  remains  constant  and  then  increases  with  the  right  well  width.  This  is  because  as 
the  right  well  width  begins  to  increase,  the  transition  energy  increases  and  the  overlap  of 
the  wave  functions  of  the  two  states  does  not  increase  very  much  and  results  in  a  decrease 
of  oscillator  strength.  As  the  right  well  width  continues  to  increase,  the  two  effects 
(increase  of  transition  energy  and  wave  functions  overlap)  are  balanced  and  oscillator 
strength  is  almost  constant.  Further  increasing  the  right  well  width,  the  increase  of  the 
wave  functions  overlap  becomes  dominant,  which  results  in  the  increase  of  oscillator 
strength. 

The  effect  of  left  well  width  on  the  energy  levels  and  oscillator  strength  is  shown 
in  figure  6.3(a)  and  (b),  respectively.  In  the  calculation,  the  right  well  and  thin  barrier 
widths  are  taken  as  70  A  and  30  A,  respectively.  It  is  seen  from  figure  6.3(a)  that  E2 
decreases  with  left  well  width,  and  Ej  is  almost  not  affected.  Compared  to  figure  6.2(a), 
the  decrease  of  E2  with  left  well  width  is  not  so  significant  as  the  decrease  of  E]  with  the 
right  well  width.  This  is  due  to  the  fact  that  the  barrier  height  for  the  left  well  is 
significantly  smaller  than  that  for  the  right  well.  The  oscillator  strength  decreases  with  the 
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left  well  width.  This  is  due  to  the  fact  that  when  the  left  well  width  increases,  the  overlap 
between  the  wave  functions  of  two  levels  decreases  and  therefore  causes  the  decrease  of 
oscillator  strength. 

Figure  6.4  show  the  effect  of  thin  barrier  thickness  on  the  energy  levels  and 
oscillator  strength.  In  the  calculation  the  right  well  width  and  the  left  well  width  are  taken 
as  50  A  and  70  A,  respectively.  It  can  be  seen  from  figure  6.4(a)  that  both  of  the  two 
energy  levels  decrease  slightly  with  the  increase  of  the  thin  barrier.  The  oscillator  strength 
also  decreases  with  the  thin  barrier  thickness,  since  increasing  the  thin  barrier  thickness 
the  wave  function  overlap  between  these  two  levels  will  decrease. 

It  is  known  from  the  single  quantum  well  case  that  more  levels  will  occur  when 
the  barrier  height  is  kept  constant  while  increasing  the  well  width.  This  also  happens  in 
this  double  quantum  well  structure.  For  example,  an  additional  level  will  appear  when  we 
increase  the  right  well  width  to  60  A  and  left  well  width  to  100  A. 

6.3.2  The  effect  of  electric  field 

As  mentioned  before,  the  transition  energy  for  this  QWIP  will  increase  with  a  blue 
shift  under  negative  bias  condition.  The  variation  of  transition  energy  with  applied  bias  is 
shown  in  figure  6.5.  As  seen  from  this  figure,  there  is  a  significant  Stark  shift  in  this 
system  with  the  applied  electric  field.  At  zero  bias  the  peak  energy  for  the  transition  is 
108.6  meV  while  at  20  kV/cm  the  peak  transition  energy  shifts  to  125.8  meV.  It  is  seen 
that  the  peak  energy  will  shift  17.2  meV  at  an  applied  electric  field  of  20  KV/cm,  which 
corresponds  to  a  peak  wavelength  shift  of  1.56  pm.  The  corresponding  voltage  drop  in 
this  double  wells  structure  is  0.15  V  at  this  electric  field,  which  is  considered  to  be  quite 
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low.  In  the  calculation,  the  right  well  width,  thin  barrier  thickness  and  the  left  well  width 
are  taken  as  50  A,  30  A  and  70  A,  respectively. 

The  wave  functions  are  less  sensitive  to  the  applied  electric  field,  as  shown  in 
figure  6.6.  As  seen  from  figure  6.6(a),  the  three  wave  functions  for  the  first  bound  level 
under  the  electric  field  of  0,  10  and  20  kV/cm  are  almost  indistinguishable.  With  the 
increase  of  applied  electric  field,  the  wave  function  of  the  second  bound  level  shifts 
gradually  to  the  right  side,  which  will  result  in  the  enhanced  overlap  of  the  wave 
functions  of  the  two  levels,  as  shown  in  figure  6.6(b).  This  is  consistent  with  the  weak 
coupling  scenario. 

As  a  result  of  the  enhanced  overlap  of  the  wave  functions  of  the  two  levels,  the 
absorption  coefficient  of  the  system  will  increase.  We  have  calculated  the  absorption 
coefficient  of  this  system  under  different  biases,  and  the  result  is  shown  in  figure  6.7.  In 
the  calculation,  the  broadening  parameter  is  taken  as  2  meV.  It  is  clearly  seen  from  this 
figure  that  the  applied  electric  field  enhances  the  absorption  in  the  system.  It  is  also  noted 
that  due  to  the  smaller  overlap  between  the  wave  functions  of  the  two  levels,  the 
absorption  coefficient  in  this  weakly  coupled  double  wells  system  is  significantly  smaller 
than  that  in  the  single  quantum  well  bound  to  bound  transition.  However,  it  is  still 
comparable  to  the  single  quantum  well  bound  to  continuum  transitions. 


6.3.3  Device  performance 

Tunneling  probability  is  an  important  quantity  in  determining  the  dark  current  of  a 
QWIP.  Figure  6.8(a)  and  (b)  show  the  tunneling  probability  of  the  right  and  left  well. 
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respectively.  It  is  seen  from  this  figure  that  with  the  increase  of  electric  field,  the 
tunneling  probability  will  increase  due  to  the  lowering  of  barrier  height.  Also,  the 
tunneling  probability  of  the  left  well  is  much  larger  than  that  of  the  right  well  due  to  the 
very  small  thickness  of  the  thin  barrier. 

The  dark  current  for  this  QWIP  structure  was  calculated  at  77  K  and  40  K,  and  the 
results  are  shown  in  figure  6.9.  In  the  calculation,  we  have  assumed  the  same  mobility 
and  saturation  velocity  for  electrons  in  both  wells,  with  //^  =  ///  =  1000  cm^Vs"'  and  = 
v^/  =  5x10®  cm/s,  and  the  device  active  area  is  210x210  pml  In  a  conventional  bound  to 
bound  transition  QWIP,  the  applied  bias  voltage  is  usually  very  high  so  that  electrons  can 
escape  efficiently  from  the  quantum  well.  As  a  result,  large  dark  current  can  be  expected 
under  high  bias  condition.  In  the  double  quantum  wells  structure,  the  second  level  is 
close  to  the  barrier  height.  It  can  be  seen  from  figure  6.8  that  even  at  a  bias  of  20  KV/cm, 
the  tunneling  probability  for  both  the  right  well  and  left  well  is  quite  high  at  E2  level 
which  means  that  the  electron  can  tunnel  out  of  quantum  wells  efficiently.  This  will  result 
in  a  reduced  dark  current  compared  to  the  conventional  bound  to  bound  transition  QWIP. 

Figure  6.10  shows  the  calculated  responsivity  of  this  QWIP.  In  the  calculation, 
the  escape  probability  is  teiken  to  be  0.5,  and  the  optical  gain  is  assumed  to  be  0.6.  The 
number  of  quantum  well  period  is  taken  to  be  6  and  /  =  4500  A.  It  can  be  seen  from  this 
figure  that  the  responsivity  for  this  QWIP  structure  is  not  very  high  due  to  the  weak 
coupling  of  the  two  wells.  However,  the  wavelength  tunability  of  this  QWIP  is  quite 
large,  which  should  be  useful  for  the  tunable  LWIR  QWIPs. 
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The  Stark  shift  effect  in  an  asymmetric  weakly  coupled  double  wells  system  has 
been  depicted  in  this  section.  It  was  shown  that  the  transition  peak  wavelength  of  this 
coupled  QWIP  structure  is  very  sensitive  to  the  applied  electric  field.  Therefore  one  can 
realize  large  Stark  shift  with  relatively  small  bias.  The  overlap  of  the  wave  functions  of 
the  two  bound  states  is  relatively  small,  but  the  absorption  coefficient  of  this  system  is 
still  comparable  to  that  of  a  bound  to  continuum  transition  in  a  single  quantum  well. 
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ENERGY  LEVELS  (eV) 


Figure  6.1  Schematic  energy  band  diagrams  of  the  asymmetric  weakly  coupled 
double  quantum  wells  at  (a)  zero  bias  and  (b)  negative  bias. 


Figure  6.2  Variation  of  energy  levels  and  oscillator  strength  associated  with  the  £7  to  £2 
transition  with  the  right  well  width.  The  left  well  width  is  set  to  be  70  A  and 
the  thin  barrier  width  is  set  to  be  30  A. 
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Figure  6.3  Effect  of  left  well  width  on  the  energy  levels  (a)  and  oscillator  strength  (b).  In 
the  calculation  the  right  well  width  is  set  to  be  50  A  and  the  thin  barrier 
thickness  is  set  to  be  30  A. 


THIN  BARRIER  WIDTH  (A)  THIN  BARRIER  WIDTH  (A) 

(b) 


Figure  6.4  Variation  of  energy  levels  and  oscillator  strength  with  the  thin  barrier 
thickness. 
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Figure  6.5  Variation  of  the  energy  levels  with  applied  electric  field.  In  the  calculation, 
the  right  well  width,  thin  barrier  thickness  and  left  well  width  are  taken  as  50 
A,  30  A  and  70A,  respectively. 


Figure  6.6  Effect  of  electric  field  on  the  wave  functions  of  ground  state  (a)  and  excited 
state(b).  In  figure  6(a),  the  ground  state  wave  functions  with  the  applied 
electric  field  of  0, 10  and  20  KV/cm  are  almost  indistinguishable.  In  figure 
6(b)  the  applied  electric  field  is  0, 10  and  20  KV/cm,  respectively  for  the 
excited  state  wave  functions  labeled  1,  2  and  3. 
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TUNNELING  PROBABILITY 


7  8  9  10  11  12  13  14 


Wav  e  le  n  g  th  (|i  m) 

Figure  6.7  Absorption  coefficient  of  the  double  quantum  well  system  under  different 
electric  fields. 
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(a)  (b) 

Figure  6.8  Tunneling  probability  for  the  right  well  (a)  and  the  left  well  (b).  In  both  (a) 
and  (b),  the  applied  electric  field  are  2, 10  and  20  KV/cm  for  the  curves 
shown  from  left  to  right. 
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Figure  6.9  Calculated  dark  current  of  the  double  quantum  well  QWIP  structure 
at  T=77  K  and  T=  40  K. 
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Figure  6.10  Calculated  responsivity  of  the  double  quantum  well  QWIP  structure 
under  different  biases. 
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